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Path entanglement is a key resource for quantum metrology. Using path-entangled states, the standard 
quantum limit can be beaten, and the Heisenberg limit can be achieved. However, the preparation and 
detection of such states scales unfavourably with the number of photons. Here we introduce sequential path 
entanglement, in which photons are distributed across distinct time bins with arbitrary separation, as a 
resource for quantum metrology. We demonstrate a scheme for converting polarization 
Greenberger-Horne-Zeilinger entanglement into sequential path entanglement. We observe the same 
enhanced phase resolution expected for conventional path entanglement, independent of the delay between 
consecutive photons. Sequential path entanglement can be prepared comparably easily from polarization 
entanglement, can be detected without using photon-number-resolving detectors, and enables novel 
applications. 

Quantum metrology, which harnesses the collective nature of entangled particles, is an emerging tech- 
nology that promises to revolutionise our way of measuring and to radically alter the limits for parameter 
estimation 14 . The most common example, which is also the one with the most technological applica- 
tions, is the estimation of the relative optical phase between two arms of an interferometer: this problem is directly 
related to the realization of atomic clocks 5 , the investigation of biological samples 6 , and the detection of gravita- 
tional waves 7 ' 8 . Generating light in an entangled state delocalised over the two arms can achieve a better scaling of 
the precision with the number of photons than is possible with any classical optical technique. Indeed, for a 
classical coherent state with average photon number N, the precision in phase estimation will mostly be limited by 
the Poissonian fluctuations of the photon number, the so-called shot noise, which contributes a scaling propor- 
tional to l/ \/N on the error of the phase estimate (this is the 'standard quantum limit'). In the quantum domain N 
photons can be prepared in an entangled superposition of all photons occupying one arm, and of all photons 
occupying the other arm. Since all photons travel together, phase is accumulated more rapidly, which allows for a 
phase estimate with enhanced precision whose error scales as 1/N, the Heisenberg limit 9 . 

The key elements for realising such an advantage are the ability to produce iV-photon path-entangled states, 
and the capability of resolving the number of photons at both output ports of the interferometer. However, scaling 
up from the few-photon level has proven to be technically challenging, since the best photon sources available 
produce multiphoton states with unwanted components, which degrades the quantum advantage 10 " 16 . This 
motivates research into less technically demanding schemes that might retain the enhanced precision of con- 
ventional path-entangled states. Here we propose breaking up iV-photon entangled states so that the entangled 
photons are distributed over AT time bins, with only one photon interacting with the sample at a time. This type of 
sequential path- entangled (SPE) state has a number of advantages. First, as we will show, SPE states can be 
produced directly from Greenberger-Horne-Zeilinger (GHZ) states 17 ; the current state of the art for producing 
GHZ states 1819 is more advanced than for NOON states 15 ' 16 , so that larger SPE states can be synthesised than NOON 
states. Second, detecting photons one at a time can be done with conventional avalanche photodetectors without 
requiring photon-number-resolving detection. Finally, SPE states sample a phase over a range of times, which 
enables the extraction of Fourier components of a dynamical phase. 

A legitimate concern with a proposal to break up and separately detect an entangled state in this way is that the 
metrological power of a collective measurement on the entangled state is lost. However it is known that strategies 
based on collective measurements offer no advantage for quantum metrology 1 , and this is why SPE states can be 
just as powerful for metrology as conventional path-entanglement. In this paper we demonstrate the key features 
of this proposal by recording the interference produced by two -photon entangled states as a time delay between 
the photons is introduced. The results show that phase super-resolution is maintained independent of the delay, 
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confirming the power of SPE states as a resource for quantum met- 
rology. We generate our SPE states by feeding a polarization-entang- 
led photon pair into a modified Mach-Zehnder interferometer 20 
(MZI). More generally, we establish that higher order SPE states 
could be deterministically produced from an initial resource of 
GHZ polarization- entangled states by means of optical switching. 

Results 

Sequential path entanglement. The archetypal path-entangled state 
is the NOON state 3 ' 21 , |M)::0N> = \N) a \0) b + ^|0) a |N) b (ignoring 
normalisation here, and below), in which AT photons are superposed 
between all propagating in one arm, a, of an MZI, and all propagating 
in the other arm b (see Fig. 1 (a)), which accrues a phase rotation N(/), 
where <j) is the relative phase between the two arms of the MZI. The 
advantage of this path-entangled state for metrology can be under- 
stood using the concept of the photonic de Broglie wavelength 22 . The 
N photons are treated as an effective single entity, and thus have an 
effective de Broglie wavelength klN i with X the single-photon 
wavelength. 

The scheme that is most commonly used to produce approximate 
NOON states is shown in Fig. 1 (a). Two Fock states \N/2) are made to 
interfere on the first beam splitter (BS) of an MZI. For N = 2, Hong- 
Ou-Mandel (HOM) interference 23 directly generates the required 
photon-bunched |20::02) state 16,24 . Unfortunately, for N > 2, gen- 
erating a genuine NOON state by HOM interference in this way no 
longer works: the resulting state is instead a superposition of \N0::0N) 
with an undesired contribution |Q(N)), which decreases the assoc- 
iated measurement precision 13,25 ' 26 . For instance, Nagata et al 
prepared a 4-photon NOON state with |Q(4)) oc |22), by interfer- 
ing 2-photon Fock states 13 . States of this form are known as 



Holland- Burnett 27 states, and, while they are interesting for lossy 
phase estimation 28 , they cannot achieve the Heisenberg limit. It has 
been proposed that |Q(N)) can be suppressed by using feedforward 29 , 
however, this proposal has not been implemented yet. 

Multiphoton SPE states can be generated without suffering from 
the complications described above that make NOON states difficult to 
produce for N > 2. We begin with a polarization- entangled GHZ 
state 1819 of N photons, 

|GHZ(iv)> = |h> 1 |h> 2 . . . \h) n +\v) 1 \v) 2 . . . \V) N , (1) 

in which |H),- (|V),-) denotes a horizontally (vertically) polarized 
photon in spatial mode i = 1, 2, N. The conversion of this 
GHZ state into an SPE state is schematically shown in Fig. 1 (b). 
Each photon from the GHZ state is delayed with respect to the next 
and then all are transferred into a single spatial mode. In our experi- 
ments we combine the photons probabilistically, using a beam split- 
ter, but in general this step could be done deterministically by means 
of a fast optical switch, such as a Pockels cell 30 . Finally, path entan- 
glement is created by directing the photons one by one onto a polar- 
izing beam splitter followed by a half- wave plate at one output mode, 
which applies the operation \H) \l) a \0) b ; |V) -» \0) a \l) b . This 
recipe allows the deterministic conversion of GHZ entanglement 
into SPE states within the MZI of the form 



|spe(n)>= nji)^ 10}^.+ njo}^ 



(2) 



where now each temporal mode tj contains a single path-entangled 
photon. Each of the AT photons in mode b picks up the interferometer 
phase and so this state also exhibits super- resolution, like the 
NOON state (see Methods). However, the time separation of the 




(|N> a |0> b +|0> a |N> b ) + |O(N)> 





nii> a( ^o> M+ n|o> a ji>^ 



time 
— ► 



Figure 1 | Preparation of NOON states and SPE states, (a) Approximate preparation of NOON states by interfering Fock states, (b) Deterministic 
preparation of SPE states converted from polarization-entangled GHZ states in a sequential way. BS: beam splitter; PBS: polarizing beam splitter; 
HWP: half-wave plate; APD: avalanche photodetectors; PNRD: photon-number-resolving detectors; (j): phase shift. 
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photons means that photon-number resolving detectors are not 
required to identify the iV-photon events containing the enhanced 
metrological signal. In principle, 'on-off detectors such as commer- 
cially-available avalanche photodetectors can be used, since there is 
only a single particle in each time bin 31 . Additionally, if the MZI 
phase is varying in time — as would be typical in pump-probe type 
experiments — the Fourier components of the phase evolution can 
be accessed efficiently via the transmission statistics of SPE states by 
tuning the time delay T = t j+ x — tj into resonance with a component 
of interest (see Methods). 

Experimental setup. Our experiment is shown schematically in 
Fig. 2. It contains two parts: the photon source (Fig. 2 (a)) which 
produces polarization- entangled photon pairs, and the modified 
MZI (Fig. 2 (b)). In the experiment, as shown in Fig. 2 (a), type-II 
spontaneous parametric down conversion 32 (SPDC) is employed to 
generate a polarization- entangled state \H) P \H) P + | V) p \ V) Pl (see 
Methods). A translation stage with um precision adds an adjustable 
delay AL X into path P l5 and then P x and P 2 are combined at a 50 : 50 
beam splitter BS1 and the resulting polarization- entangled state is 
coupled into an optical fibre in mode P 4 for delivery to the MZI. We 
note that BS1 passively combines the photons into the same spatial 
mode in a probabilistic way, which allows us to study several delay 
settings, including AL X = 0. But for a sufficiently large delay AL X 
where £ is the single photon coherence length, the photons are well- 
separated and could be deterministically combined using a fast 
optical switch, allowing the efficient conversion of higher-photon- 
number GHZ states into SPE states. The state emerging from the 
fibre is converted into an SPE state of the form (2) at a polarizing 
beam splitter, and a nm-precision piezoelectric actuator mounted on 




a urn-precision translation stage is used to control the optical path 
length difference AL 2 within the modified MZI. 

To find the position AL 2 = 0 we first maximise the visibility of 
single-photon interference in the MZI (see Fig. 2 (d)). For these 
measurements BS1 is replaced by a polarizer oriented at 45°, so that 
photons emitted from the source into path P x are delocalised at the 
polarizing beam splitter, producing a fringe pattern in the singles 
count rates recorded in the paths P 7 , P 8 , as AL 2 is scanned (photons in 
path P 2 are simply dumped). The full width at half maximum of the 
pattern provides a measurement of the single-photon coherence 
length £ ~ 130 /mi, consistent with the result deduced from the 
measured HOM interference dip shown in Fig. 2 (c) (see 
Methods). We then set AL 2 ~ 0 and scan the piezoelectric actuator 
in steps of 10 nm to observe the fringe spacing, which matches the 
single photon de Broglie wavelength (simply the optical wavelength) 
X « 810 nm (see Fig. 3 (a)). 

Next, we verify phase super-resolution with a standard NOON state 
by setting AL X = 0, so that both photons emitted from the source are 
temporally overlapped in the MZI. As shown in Fig. 2 (e), the coin- 
cidence count rate at detectors D x & D 2 exhibits an interference 
pattern with a broader envelope than the single photon interference 
pattern. This is due to spectral entanglement of the photon pairs 
emitted by our source, which introduces temporal correlations that 
are limited only by the coherence time of the downconversion pump 
laser, extending significantly beyond the single-photon coherence 
time (see Methods for details). We comment that to our knowledge 
this is the first time the extended interference pattern of a |20::02) 
state generated with a spectrally anti- correlated photon pair source 
has been fully mapped out 24 . The fine-grained scan of the region 
around AL 2 = 0 in Fig. 3 (b) shows the expected half-period fringes 



12000 




AL 2 (um) 

Figure 2 | Schematic experimental setup and one- and two-photon interference, (a) Preparation of coaxial polarization- entangled state. Polarization 
entangled photon pairs are generated by type II degenerate SPDC in a nonlinear crystal. BS 1 serves as a probabilistic optical mixer coupling both entangled 
photons into path P 4 . AL X is tuned by scanning the position of a prism mounted on a translation stage, (b) A modified MZI. (c) The two-photon 
coincidence rate measured in paths P 3 and P 4 exhibits a characteristic dip as AL X is scanned, indicative of HOM interference at BS1. (d), (e) One- and two- 
photon interference patterns. BBO: /^barium-borate crystal; COMP: HWP and BBO (1 mm); POL: polarization plate; PIEZO: piezoelectric actuator. 



SCIENTIFIC REPORTS | 3 : 1779 | DOI: 1 0.1 038/srep01 779 



3 




0.8 1.0 1.2 

AL 2 (|jm) 

Figure 3 | Characteristics of one- and two-photon interference, (a) Single-photon count rate in path P 7 (red solid circles) and at P 8 (blue solid circles) as 
the delay AL 2 is varied. Both output modes of the interferometer are directed to similar detection set-ups, which provides a consistency check, (b) Two- 
photon coincidence rates at detectors D u D 2 (red solid circles) and at detectors D 3 , D 4 (blue solid circles) versus AL 2 ; events with a time separation within 
ALJc are considered as coincident. 



associated with the reduced de Broglie wavelength of the |20::02) 
state. 

Independence of the separation between photons. To verify that 
SPE states have the same reduced de Broglie wavelength as the NOON 
state, independent of the delay between successive photons entering 
the MZI, we carried out experiments for AL X = 200 fim and 
10,000 |im. The setting AL X = 200 fjm > £ is larger than the 
single-photon coherence length but still less than the coherence 
length of the pump laser £ pump ~ 300 fim, while AL 1 = 
10,000 /nm is significantly larger than £ pump . Fig. 4 shows the 



interference patterns in the coincidence rates for detectors D Y and 
D 2 produced by scanning AL 2 . For each of the scans, the periodicity is 
always about 405 nm, half of the single-photon wavelength X = 
810 nm, which confirms that SPE states retain the phase super- 
resolution of a conventional NOON state, even when consecutive 
photons are separated well beyond their individual coherence 
lengths. It is clear that the coincidence rates fall as AL X is 
increased, but this simply represents a reduction in the SPE prepa- 
ration efficiency. This is because for AL X < HOM interference 
at BS1 boosts the probability that photon pairs are coupled into 
our MZI, whereas for ALi>150^m there is no longer HOM 



400 



400. 
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Figure 4 | Two-photon interference envelope and fringes versus interval AL 2 . Error bars are given by Poissonian statistics, (a), (b) AL X = 0 urn. 
(c), (d) ALi = 200 ^m. (e), (f) AL Y = 10, 000 ^ni. 
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interference. In addition, the fibre coupling efficiency drops for large 
AL j. In particular, the setting ALj = 10,000 exceeds the Rayleigh 
range defined by our fibre, causing a significant drop in the 
coincidence rate for this measurement. However, the peak 
visibilities of the patterns in parts (b)-(f) of Fig. 4 remain high 
throughout the scan, at 0.98 ±0.01, 0.95 ± 0.03 and 0.98 ± 0.03, 
respectively. From this we conclude that the quality of the bi-photon 
interference is independent of AL l5 and the reduced photonic de 
Broglie wavelength is independent of the distance between photons 
in an SPE state. 

Discussion 

We have introduced SPE states as an alternative resource for 
quantum-enhanced optical interferometry. Conventional NOON 
states are difficult to produce, and photon-number-resolving detec- 
tors 33 " 35 are required to extract a metrological signal 36 . While multi- 
plexed on-off detectors can be realized easily, they only count photon 
number approximately, and are therefore not suitable for metrolo- 
gical applications. By contrast, SPE states can be generated determi- 
nistically from GHZ states, which are easier to produce than NOON 
states. Photon-number-resolving detectors are also not needed, since 
photons can be detected sequentially 1 . In practice, the performance 
of photodetectors is critical to achieving a metrological advantage 
with NOON states 28 , and this remains the case for SPE states also. 
However, realistic quantum interferometry could be realised using 
recently developed highly efficient, fast, detectors 37 . 

Many experiments have shown that two photons arriving simul- 
taneously make a substantial difference in second harmonic genera- 
tion 38 , Kerr-like interactions 39 and Rydberg atomic gases 40 . This 
detrimental effect makes phase estimation modeling complicated. 
Furthermore, quantum metrology is often motivated by the need 
to probe refractive index variations in samples with extreme sens- 
itivity to optical radiation 41 . NOON states can excite 2 -photon or 
multiphoton transitions in a sample that scatter light out of the 
interferometer, and such transitions can ionise quantum systems 
and damage biological molecules. The influence of multiple photons 
will become stronger if more photons are present at one time, as 
could be the case with high-iV NOON states, which are the ultimate 
goal of this avenue of research. By contrast even large-iV SPE states 
can excite only single-photon transitions because the constituent 
photons arrive one at a time. 

It is also notable that SPE states can be used to extract the har- 
monic spectrum of a time-varying phase: if the consecutive time bins 
in an SPE state are separated by a constant delay T, the phase is 
sampled periodically, producing a signal proportional to the ampli- 
tude of the Fourier component of the dynamical phase at a frequency 
/= l/T (see Methods). 

We experimentally confirmed that 2 -photon SPE states posses the 
same phase super-resolution as 2-N00N states. Unlike schemes for 
achieving super- resolution via post-selection on classical light 42 , this 
super- resolution is of quantum origin, and SPE states therefore rep- 
resent a convenient and powerful platform for quantum-enhanced 
phase sensing via optical interferometry. 

Methods 

Preparation of polarization-entangled photonic source. As shown in Fig. 2 a CW 
semiconductor laser produces a blue beam (power 34.5 mW, waist 100 /mi, central 
wavelength 405 nm) that pumps type-II degenerate SPDC in a 2 mm jS-barium- 
borate (BBO) crystal, producing pairs of photons with central wavelength 810 nm 32 . 
The down- converted signal and idler photons have different propagation velocities 
and travel along different paths inside the crystal due to its birefringence. The 
resulting walk-off effects are compensated by a combination of a half- wave plate and a 
1 mm BBO crystal in each arm. Photons are then coupled into single mode fibres and 
we prepare a polarization-entangled state \H) P \ V) Pi + 1 V) p \H) Pi . The state is 
further transformed to be \H) P \H) P + 1 V) p \ V) Pi by placing a half- wave plate with 
45° in mode P 2 . The singles count rates measured directly in each of P l5 P 2 are about 
80,000 per second, and the coincidence rate is about 20,000 pairs per second. The 
visibilities for the polarization correlations are about 98% in the \H)/\V) basis and 93% 



for the | + 45°)/ 1 — 45°) basis. This shows that we have produced a high-quality 
polarization-entangled photonic source. 

Success rate for combination of photon pairs at BS1. When AL X > £, the photons 
do not overlap at BS1 and each photon has the same probability to be transmitted or 
reflected. Thus, the photon state in P 4 can be a single-photon state |H> or \ V) with 
probability 0.25, respectively, or the maximally polarization-entangled state with 
probability 0.25 (the rest of the time, P 4 contains no photons). When AL X < HOM 
interference suppresses the single photon probabilities, boosting the probability for 
coupling both photons into P 4 up to 1/2 for AL X = 0. HOM interference with a 
visibility of 94.5 ± 0.4% is shown in Fig. 2 (c), obtained by scanning AL X and 
measuring the coincidence rate at two output modes P 3 and P 4 (Fig. 2 (a)). The high 
visibility certifies good spatial and temporal overlap for the entangled photon pairs 
emitted by our source. A Gaussian fit to the HOM dip in Fig. 2 (c) yields the single- 
photon coherence length £ ~ 126 fim. 

Single- and bi-photon coherence time. Here we briefly outline the origin of the 
differing widths of the single and bi-photon interference patterns shown in parts 
Fig. 2 (d) and (e). We consider that our downconversion source creates photon pairs 
in the entangled state 

\\J/) = ^£(co,co')\co + coq,co' + coo) dcodco', 

where £ is the joint spectral amplitude for emitting signal and idler photons with 
frequencies co + co 0 , co' + co 0 , with co 0 the centre frequency of the downconversion. £ 
generally extends over a broad bandwidth set by the phasematching conditions in the 
downconversion crystal, but energy conservation enforces tight frequency anti- 
correlations, with a correlation width equal to the pump laser bandwidth. 

For the single photon interference measurement (Fig. 2 (d)) we dump photons P 2 , 
leaving the photons in P 1 in the mixed state p = §f(coi,co 2 )\coi +co 0 ) 
(co 2 -\-co 0 \dco 2 dco 2 , where /( co i,co 2 ) = J C(coi,co)[*(co 2 ,co)dco, The MZI implements 
the transformation 

P ^ P>= J/( (yi ' W2 )[ 1±ei(ftJl+aJ ° )AL2/C ] [l±e- i(fU2+fUo)AL2 / c ] 

I coi + cqq ) (co 2 + a» 0 1 do)\ dco 2 

at the two output ports P 7 , P 8 (postselected on no photon in the other port), and the 
single photon count probability is P = tr{ p' J" dft>|co) (co\ } = 

Jp(w)|l±e i ( ffl+ffl °) Al2 / c | dco, where p{co) = j{co, co) is the marginal spectrum. This 
exhibits fringes with period X = Incojc and visibility V given by a Fourier transform, 
V(AL 2 ) = J r co {p(co)}(AL 2 /c), so that the coherence length £ is determined by the 
inverse of the bandwidth of p(co). Note that the fringe patterns for the two output 
ports are n out of phase, as seen in Fig. 3 (a). 

For the measurement of the two photon NOON state interference (Fig. 2 (e)), the 
MZI produces the state 

\xj / ') = |cKwO[l + e i(ffl+w ' +2wo)Ai2/c ]|co + wo,w / + wo>dwdw / 

at each of its two output ports (postselected on no photons in the other port). The 
coincidence probability at detectors D u D 2 is then 

V = tr { | {if/' |J dco dco' | co,co') (co,co r | }, which produces an interference pattern with 
period XI2 and visibility given by V(AL 2 ) = g(AL 2 /c, AL 2 /c), where 
g(t,t') = J- a),a)' { \£(co,co') | 2 } (t,f) is the two-dimensional Fourier transform of the 
joint spectral intensity \C(co, co')\ 2 . The tight frequency anti-correlations between the 
photons produce strong temporal correlations such that the width of g(AL 2 /c, AL 2 /c) 
is proportional to the coherence length £ pump of the pump laser. 

Super-resolution with SPE states. To convert NOON states into SPE states we simply 
introduce a delay AL X between photons, which for the bi-photon case corresponds to 
the transformation £(co,co') -»£(co,(o')e i ( co + < °o)ALi/c < However, as described above, the 
bi-photon interference visibility depends only on the joint spectral intensity, which is 
unaffected by this phase rotation. This is why our SPE state retains phase super- 
resolution for all values of AL X . For higher photon number states with N > 2 the 
visibility is determined by higher order correlation functions, but the coincidence rate 
is always invariant under phase rotations of the joint spectrum, so that the 
metrological power of SPE states is independent of the delay between photons. 

Extracting dynamical phase information with SPE metrology. An interesting 
feature of SPE states is that the photons sample the MZI phase over a set of discrete 
times tj y meaning that the metrological signal is of the from V cos 2 (0), where Vis a 
visibility, and where O = Ylf= Y (j) (tj) . Substituting in the Fourier relation 
(j)(t) = (2n)~ 1 ^ 2 ^(co)e 10Jt dco, where 0 is the spectrum associated with the time- 
varying phase (note that here co is the Fourier conjugate of t, and not an optical 
frequency), we obtain 

<P = [AH^Mcko; AH = (2ft)-'/W*->^ Sin( f» r / 2) , (3) 
J sm(coT/2) 

where we assumed a constant inter-photon delay tj = t 0 + (j — 1 ) T (in our experiment 
we had T = ALJc). The function A comprises a series of spectral peaks of height 
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(2n) l/2 e i(0kto , located at the frequencies co k = InkIT, k = . . ., -2, - 1, 0, 1, 2, each 
with width Sco = 4n/NT. In the limit that AT is large, we may approximate A as a series 
of Dirac delta functions, so that the SPE signal can be written as 

co 

Ooc e~ i( ° kk h- (4) 

k= — co 

We would like to extract information about the Fourier components fa = (f>(a)k) of 
the phase; as an example we consider the quantity D T = <&\ to =0 — <D| fe = T , where t = 
nlwi = 772, calculated from two measurements with different values of the initial 
arrival time t 0 . The zero -frequency component fa is removed from D T , and we have 

D T (x5Re{^}, (5) 

where we have used the fact that ^>_ k = fa (since 0 is real), and where we also assumed 

that the phase variation is sufficiently narrowband that | fa 1 <C | fa | , which will hold for 

a wide class of smooth phase dynamics with compact temporal support of duration 
T^>T. This shows that just two measurements with SPE states provide access to the 
Fourier component of a dynamical phase at a frequency/ = co^ln = II T fixed by the 
delay T between the probe photons. We note that the accuracy of this procedure is 
ultimately limited by the timing jitter of the nonclassical light sources that produce 
the photons in the SPE state, but parametric photon sources have demonstrated 
extremely low timing jitter 43 . 
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